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Abstract 

We investigate the polarization dependence of optical nonlinearity enhance- 
ment for a uniaxial anisotropic composite of metal nanocrystals in a dielectric 
host. Three cases are distinguished depending on whether the polarization is 
parallel, perpendicular or unpolarized with respect to the axis of anisotropy. 
For the parallel polarization, the results show that the 3D results are qualita- 
tively similar to the 2D case reported recently. For the perpendicular polar- 
ization, the results are markedly different from the parallel counterpart: In 
contrast to the absorption, the enhancement factor actually increases with the 
anisotropy. Thus the separation of the absorption and enhancement peaks be- 
comes even more pronounced than the parallel polarization case. These results 
indicate a strong polarization dependence of the nonlinear optical response. 
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I. INTRODUCTION 



With the advent of high-power coherent light sources, optical nonlinearity has become a 
rapidly growing field both in basic research as well as in practical applications [ffj. The use- 
fulness of optical nonlinearity depends on the availability of materials with a large nonlinear 
susceptibility and a short response time. While the search for a homogeneous material that 
possesses a large nonlinear susceptibility continues, the use of composite materials has been 
proposed as an alternative for enhancing the optical nonlinearity, through the inhomogene- 
ity of the local electric field and/or the surface-plasmon resonances These enhancement 
mechanisms are sensitive to the composite microstructure. 

In a recent work , we proposed the use of structurally anisotropic composites to enhance 
the optical nonlinearity. Anisotropy is clearly a phenomenon common in most materials, 
and in optical phenomena. Anisotropy can be an intrinsic material property or can be 
induced, e.g., by the application of fields. The nonlinearity can be enhanced by using 
nonspherical (e.g., ellipsoidal) particles 0. Induced nonlinearity such as the electrostrictive 
mechanism (through gradients in EM fields) that operates in the microsecond regime was 
recently reviewed [|J. 

In Ref. 0, we considered composites made of components which are locally isotropic in 
their material properties, but with electric-field induced anisotropic microstructures. Here 
the characteristic length scale of the microstructure is assumed to be much less than the 
wavelength, so that the optical response of the system is characterized by an effective dielec- 
tric tensor. We showed that by controlling the magnitude of the applied electric field, for 
example, one can vary the degree of geometric anisotropy and thereby tune the system to 
achieve maximum optical nonlinearity in accordance with the properties and volume frac- 
tion of nonlinear materials 0. Furthermore, we proposed the use of the electrorheological 
(ER) effect |||| to realize the desired geometric anisotropy, e.g., during the annealing pro- 
cess, in material systems containing a nonlinear optical component. An important benefit 
of the anisotropic microstructures considered in this work and Ref. is the separation of 
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the nonlinear enhancement peak from the absorption peak, thus raising the application 
potential of the microstructure enhancement mechanism. 

While preliminary results were presented for the parallel polarization in 2D ||, it is 
also necessary to obtain the results for the perpendicular polarizations and hence results 
for the unpolarized response. Moreover, it is also instructive to examine the more realistic 
3D case. In what follows, we shall use the spectral representation of local fields |7J to 
model the anisotropy enhancement effect. The spectral representation offers the advantage 
of separation of material parameters from the microstructure information, contained in the 
spectral density function m(s). It has been shown recently that the optical nonlinearity 
enhancement effect is directly linked to the behavior of m(s) fllO |. 

The paper is organized as follows. In the next section, we present the spectral represen- 
tation of anisotropic composite media. In section III, we present results for the absorption 
and enhancement peaks in the optical response, followed by a summary of our results. 



II. SPECTRAL REPRESENTATION OF ANISOTROPIC COMPOSITE MEDIA 

We consider a composite metal-dielectric system, with metal being the nonlinear optical 
component. We assume the particle size is much smaller than the wavelength of light so 
that the quasi-static approximation can be used. The local constitutive relation is given 
by D = (e + x|E| 2 )E, where e is the (position dependent) dielectric constant and \ is 
the (position dependent) third-order Kerr nonlinear susceptibility. We model this material 
system by a lattice model, with assigned probability of metal particle occupation at each 
lattice site. To introduce uniaxial anisotropy, we assign two different occupation probabilities 
parallel (pn) and perpendicular (p±) to a particular direction of an occupied site. We choose 
several values of the anisotropy ratio r = p±/p\\, varying from isotropic (r = 1) to highly 
anisotropic (r = 0). We invoke the effective- medium approximation (EMA) for anisotropic 
media The coupled EMA self-consistency equations read 
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where z\\ and z± are the parameters parallel and perpendicular to the applied dc field. In 
Eqs.(jl]) and (fj), e% denotes the dielectric constant of metal, and e 2 that of the dielectric 
component; e\\ and e± denote, respectively, the parallel and perpendicular components of 
the effective dielectric tensor. In 2D, the ^-parameters are given by M 



tan 1 \/e_i_/ei| 
tan Je\\/e± 



tan 1 \/e\\/e± 
tan Je±/eu 



Let U = ej_/ en . In 3D, the ^-parameters are given by 



tan' 1 y/2U + U 2 

tan" 1 V2U + U 2 /U 

z± = ; / 6 

tan" 1 U/V2U + U 2 V ; 

The above theory is applicable to symmetric microstructure as in granular composites, and 
is an extension of the Bruggeman theory to anisotropic composite media |§. The depolar- 
ization factors are related to the ^-parameters through D\\ = l/(l + z\\) and D± = l/(l + z±). 
There is a simple relation for the depolarization factors for uniaxial anisotropy in d dimen- 
sions 0: 

D\\ + (d-l)D ± = 1, 

which can be readily verified by using Eqs.(^)-(Bp. For the dispersion microstructure, as in 

colloidal suspensions, it is more appropriate to use the Maxwell-Garnett approximation |§: 

ey - e 2 e x - e 2 



e\\+z\\e 2 ei + z\\e 2 ' 



(7) 



— =P^—_ • ( 8 ) 

Below we consider only the EMA case. For two-component composites, it has proved 
convenient to adopt the spectral representation of the effective linear response 0: Let 
v — 1 — ei/e 2 , w a = 1 — e a /e 2 (a = \\ or _L), and s = l/v, we find 




where m a (s') is the spectral density which is obtained through a limiting process: 

n^a(s') = lim Im w a (s' + irf). (10) 

In what follows, we consider the 3D case. The EMA can readily be solved in the spectral 
representation. In our numerical calculations, we choose the real part at several hundred 
equally spaced values across the interval < s' < 1, and the imaginary part rj to be some 
small positive value. The actual value of rj is unimportant. We found that rj = 0.001 gives 
acceptable results by checking the sum rule: 

f m a (s')ds' =p a . (11) 

JO 

From any given value of s, the ratio ei/e 2 is calculated from s = s'+irj, and subsequently e a , 
hence w a , may be evaluated numerically. In what follows, we denote p» = p for convenience. 

For isotropic composites, r = 1, z\\ = z± = 2 and the EMA self-consistency equations 
can be solved analytically. The spectral density is given by JTIJ : 

m(s',p) = ~(3p- l)e(3p- l)5{s') + JL^/( S '- S1 )( S2 - S '), (12) 

valid for si < s' < S2, where 

si = ^(l+p-2^/2p(l-p)), s 2 =^(l+p + 2pp{l-p)y (13) 

We first consider the parallel response, ey. In Fig.|l], we plot the spectral density rn\\{s) 
(associated with en) against s for several values of p and r. For small p, e.g., p = 0.1, the 
metal nanoparticles form essentially isolated clusters and the spectral function exhibits a 
single peak centered around s = 0.2. More interestingly, as the anisotropy increases, the 
spectral peak becomes narrower. This can be understood by the fact that when r — > 0, more 
chain-like isolated clusters are formed, leading to a narrow spectral peak. 

When p gradually increases, the interaction among clusters occur, leading to a broad 
spectral peak in the isotropic case. However, the spectral peak remains relatively narrow in 
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the anisotropic cases. When p increases further, the metal clusters begin to percolate the 
system in the isotropic case, leading to a delta-function contribution to the spectral density 
at s — 0. As p > 1/3, we observe the familiar edge singularity near s = for the isotropic 
case. As is evident from the figures, the 3D results are qualitatively similar to the 2D ones 



spectral function substantially. This is the source of the pronounced effect on nonlinearity 
enhancement. 

We next present the perpendicular response, e±. In Fig.||, we plot the spectral density 
m±(s) (associated with e±) against s for several values of p and r as in Fig.l. For small p, 
e.g., p = 0.1, again the metal nanoparticles form essentially isolated clusters. The spectral 
density is generally smaller in magnitude. More interestingly, as the anisotropy increases 
(r — > 0), mj_(s) decreases with r. This is due to the fact that chain-like clusters are isolated 
from one another in the perpendicular direction. As r decreases, the spectral function 
splits to a two-peak structure. For larger p, the decrease of mj_(s) with r persists. Since 
the absorption is directly proportional to mj_(s), we obtain an interesting result that the 
absorption actually decreases with the strength of anisotropy. 

III. POLARIZATION DEPENDENCES OF THE ABSORPTION AND THE x (3) 

ENHANCEMENT PEAK 

When the spectral density m a (s) is known as a function of s and p, the effective linear 
response can be calculated from Eq.@. We adopt the Drude model for the dielectric function 
of metal nanoparticles: 



where u p is the plasma frequency and 7 is the damping constant. We choose 7 = 0.0lu p and 
€2 = 1.77 for our model calculation. In this work, we assume only the metallic component 
to be nonlinear. 



except for some fine details. From the results, it is clear that anisotropy can change the 




(14) 



6 



In Fig|| the absorption peak Im(e||) is plotted against frequency uj for various r and 
p. In Fig.[|, the absorption peak Im(ejJ is plotted against frequency uj for various r and p. 
As evident from the results, the absorption peak exhibits similar behavior as the spectral 
density does. This is attributed to the fact that the absorption is related to the imaginary 
part of the effective dielectric function. From the results, it is also clear that anisotropy can 
have an important effect on absorption. We next study the effect of anisotropy on the 
enhancement factor. 

If a plane-polarized electromagnetic wave of amplitude Eq with the polarization parallel 
the uniaxial anisotropy axis is incident upon the composite system, the local field averages 
are given by [[L(J 



P<l^| 2 >= C'^-T^rEl (Hi) 



10 \s — s' 

From the average local fields, we calculate the effective nonlinear response as: 



X\\\E Q \ 2 E*=p X x(\E 1 \ 2 ){Ef). (17) 



This expression results from the mean-field approximation [ID . 

In Fig.[5], we plot the enhancement factor, |xii|/xi, against uj for various r and p. Here 
we observe that anisotropy has indeed a pronounced effect on the enhancement peak, as 
expected. As p increases, the enhancement peak exhibits a red shift, in analogy to the 
Maxwell-Garnett microstructure. These results are qualitatively similar to the 2D ones 
reported recently [Qj. 

Similar considerations apply to the perpendicular polarization. In Fig.^, we plot the 
enhancement factor, |x_l|/xi, against uj for various r and p. Here, in contrast to the absorp- 
tion (Fig.|j), the enhancement factor actually increases with the anisotropy as r decreases 
from unity. For small p, the enhancement peak becomes extremely narrow at r = 0.01. 
However, the enhancement peak must vanish at r = (because p± = in this limit). The 
separation of the absorption and enhancement peaks is thus even more pronounced for the 



perpendicular polarization. Here we observe that anisotropy has indeed a pronounced effect 
on the enhancement peak, as expected. We conclude that anisotropy can enhance the non- 
linearity and its figure of merit, due to the separation of the absorption peak from the x^ 
enhancement peak. 

We can calculate the optical response of unpolarized light, which is defined as the average 
of the results of parallel and perpendicular polarizations. 

ee = ^(e||+ej_), (18) 

Xe = \{X\\+XA.)- (19) 

The factor 1/2 is due to the fact that unpolarized light has equal components in the parallel 
and perpendicular directions, which is valid for normal incidence and when the anisotropic 
uniaxial axis is parallel to the flat surface. For oblique incidence, we must deal with the s 
and p polarizations separately and compute the average. For the s polarization, i.e., when 
the electric field is normal to the plane of incidence, we expect to obtain the same results as 
those of normal incidence. For the p polarization, i.e., when the electric field is parallel to 
the plane of incidence, the results will depend on the angle of incidence. 

In Fig.[7|, the absorption peak Im(e e ) for unpolarized light is plotted against frequency 
oj for various values of r and p. In Fig|| the enhancement peak |% e |/Xi is plotted against 
frequency u for several values of r and p. It should be remarked that the parallel case does 
not contribute to the average at r = because p± = in this case. These results may be 



observed in annealed Au:SiC>2 composites [I I 



In order to show the enhancement effect more clearly, we plot in Fig|| the figure of merit 
(FOM), defined by dividing Xe/xi °f Fig-D by Im(e e ) of Fig.[7] for the unpolarized case. As 
is evident from the figures, large FOM occurs at large u. However, if we are interested in 
smaller u>, say, u/uj p < 0.5 because of the quasi-static condition, the largest FOM of about 
120 xi ( esu cm ) is achieved at p = 0.1 and r = 0. 
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DISCUSSION AND CONCLUSION 



We have investigated the polarization dependence of optical nonlinearity enhancement 
in an anisotropic composite of metal nanocrystals in a dielectric host. For the perpendicular 
polarization, the results are found to be markedly different from the parallel counterpart. In 
contrast to the absorption, the enhancement factor actually increases with anisotropy. The 
separation of the absorption and enhancement peaks is thus even more pronounced for the 
perpendicular polarization. 

We have computed the spectral density in the anisotropic effective-medium approxima- 
tion. The results support our proposal that there can be very large anisotropy-induced 
enhancement of nonlinearity. It is further proposed that such enhancement effect may be 
realized experimentally through the application of the electrorheological effect (e.g., during 
the annealing process), with the possibility of achieving even larger optical nonlinearity than 
that reported in . 



It is relevant to note that in this context that the effects discussed in the present paper 
are qualitatively different from those of a dispersion of aligned anisotropic particles (such as 
sticks) with nonlinear optical properties. Whereas in the case of aligned sticks the enhance- 
ment peak always coincides with the absorption peak, here the separation of the two makes 
the figure of merit much more attractive. Physically, the difference arises from the different 
microstructures. In the present case, the connectedness of the anisotropic microstructure is 
crucial for the predicted effects. 

Colloidal systems with nonlinear optical particles in the range of 0.1 fim are good can- 
didates for induced-anisotropy enhancement. Here the particle size range is determined by 
considerations that the particles should be smaller than the optical wavelength, yet large 
enough so that the Brownian motion is not strong enough to overwhelm the field effect (ER 
effect). 

It should be noted that when an intense dc field and an EM field are applied simultane- 
ously, there can be an enhancement in nonlinearity due to electrostriction H . However, the 
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results of the present work implies that the induced anisotropy in the microstructure can 
have an even more significant enhancement in nonlinearity. 
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FIGURES 

FIG. 1. Spectral density m»(s) of parallel polarization of the anisotropic EMA plotted as a 
function of s, for several values of r = 0, 0.25, 0.5, 0.75 and 1. (a) p = 0.1, (b) p = 0.3, (c) p = 0.5 
and (d) p = 0.9. 

FIG. 2. Spectral density m±(s) of perpendicular polarization of the anisotropic EMA plotted 
as a function of s, for several values of r = 0.01, 0.25, 0.5, 0.75 and 1. (a) p = 0.1, (b) p = 0.3, (c) 
p = 0.5 and (d) p = 0.9. 

FIG. 3. The absorption peak Im(eii) plotted against frequency uj, for several values of r = 0, 
0.25, 0.5, 0.75 and 1. (a) p = 0.1, (b) p = 0.3, (c) p = 0.5 and (d) p = 0.9. 

FIG. 4. The absorption peak Im(ej_) plotted against frequency ui, for several values of r = 0.01, 
0.25, 0.5, 0.75 and 1. (a) p = 0.1, (b) p = 0.3, (c) p = 0.5 and (d) p = 0.9. 

FIG. 5. The enhancement peak |x|||/xi plotted against frequency uj, for several values of r = 0, 
0.25, 0.5, 0.75 and 1. (a) p = 0.1, (b) p = 0.3, (c) p = 0.5 and (d) p = 0.9. 

FIG. 6. The enhancement peak |x_L l/Xi plotted against frequency uj, for several values of 
r = 0.01, 0.25, 0.5, 0.75 and 1. (a) p = 0.1, (b) p = 0.3, (c) p = 0.5 and (d) p = 0.9. 

FIG. 7. For unpolarized light, the absorption peak Im(e e ) plotted against frequency uj, for 
several values of r = 0, 0.25, 0.5, 0.75 and 1. (a) p = 0.1, (b) p = 0.3, (c) p = 0.5 and (d) p = 0.9. 

FIG. 8. For unpolarized light, the enhancement peak |xe|/xi plotted against frequency uj, for 
several values of r = 0, 0.25, 0.5, 0.75 and 1. (a) p = 0.1, (b) p = 0.3, (c) p = 0.5 and (d) p = 0.9. 

FIG. 9. For unpolarized light, the figure of merit plotted against frequency uj, for several values 
of r = 0, 0.25, 0.5, 0.75 and 1. (a) p = 0.1, (b) p = 0.3, (c) p = 0.5 and (d) p = 0.9. 
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